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We study the effect of the interfacial transparency on the Josephson current in a diffusive ferromagnetic con- 
tact between two superconductors. In contrast to the cases of the fully transparent and the low-transparency in- 
terfaces, the current-phase relation is shown to be nonsinusoidal for a finite transparency. It is demonstrated that 
even for the nearly fully transparent interfaces the small corrections due to weak interfacial disorders contribute 
a small second-harmonic component in the current-phase relation. For a certain thicknesses of the ferromagnetic 
contact and the exchange field this can lead to a tiny minimum supercurrent at the crossover between and n 
states of the junction. Our theory has a satisfactory agreement with the recent experiments in which a finite 
supercurrent was observed at the transition temperature. We further explain the possibility for observation of a 
large residual supercurrent if the interfaces have an intermediate transparency. 
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I. INTRODUCTION 



Ferromagnet-superconductor hybrid structures exhibit 
novel and interesting phenomena which have been studied 
extensively in the recent years (for a review see Refs. 1 and 
|2)- These systems provide the possibility for a controlled 
study of coexistence and competition of the ferromagnetism 
and the superconductivity. One of the most interesting effects 
is the possibility of forming the so-called tt Josephson junc- 
tion in superconductor-ferromagnet-superconductor (SFS) 
structure9^ii^i^ii2i2ii£iiLi^. In a tt junction the ground-state 
phase difference between two coupled superconductors is tt 
instead of as in the ordinary junctions'. The existence 
of the TT junction in layered SFS systems, which was first 
predicted by Bulaevskii et al.,^ and Buzdin et al^ to occur 
for certain thicknesses and exchange field energies of the 
F layer, has been confirmed in experiment"^. Ryazanov et 
al!^ observed the formation of the tt junction in diffusive 
SFS junctions by measuring nonmonotonic variations of the 
Josephson critical current on changing the temperature. At the 
transition from to tt coupling the critical cuiTent undergoes a 
sign change and thus passes through a zero-current minimum. 
This appears as a cusp in the temperature dependence of the 
absolute values of the critical current which is measured in the 
experiment. Later the nonmonotonic temperature dependence 
of the critical current in diffusive SFS contacts was observed 
in other experimentsi^tiiii^ii^ and was attributed to the — tt 
transition induced by the feiTomagnetic exchange field. The 
— TT transition have been studied theoretically by several 



authors in cleans* 



8.10.11.12 



and diffusivi 



,5.8.12.17 



F Josephson 



contacts for different conditions and strength of disorder 
at the FS interfaces. The theory turned out to be in good 
agreement with experiment!^ 

Very recently a new experiment revealed another charac- 
teristic of the — TT transition in SFS junctions, which was 
not detected before. While in the early experiments^ii^iiiii^ii^ 
the critical current at the crossover temperature was found to 
vanish, Sellier et al.^"^ reported the existence of a finite small 
supercurrent at the transition temperature To^r in diffusive SFS 
junctions for a certain thickness of the F layer. A finite critical 



current at the crossover was predicted theoretically'"" in bal- 
listic F contacts, and was attributed to the appearance of the 
second harmonic in the Josephson current-phase relation. In 
Ref- Il9. further analysis of the residual supercurrent by apply- 
ing high-frequency excitations showed half-integer Shapiro 
steps, confirming the appearance of the second harmonic in 
the Josephson cuiTent-phase relation. This observation was 
confirmed by other experiments in Ref. 20, although it was 
argued that the appearance of the half-integer Shapiro steps 
in the current-voltage characteristic of the SFS contacts may 
also originate from a non uniform thickness of the F layer. 

For diffusive SFS contacts the theoretical studies mostly 
concentrated on either the perfectly transparent FS interface 
or the low-transparency case. For these two limiting cases 
and in the weak-proximity-effect regime using the linearized 
Usadel equations^' implemented with the Kupriyanov and 
Lukichev^^ (KL) boundary conditions at the FS interfaces, the 
current-phase relation (CPR) is found to be always sinusoidal. 
Melin2^ was the first who proposed a theoretical description 
of the sin(2(p) term of the CPR in SFS junctions in which 
a strong decoherence in the magnetic alloy can explain the 
magnitude of the residual supercuiTent at the — tt transition. 
Latter Houzet et al}'^ and Buzdin^^ independently presented 
an explanation of the nonsinusoidal CPR by considering the 
effect of magnetic impurities on the Josephson current in dif- 
fusive SFS contacts using the KL boundary condition in the 
highly transparent limit. They showed that the isotropic^'' or 
the uniaxial^^ magnetic scattering in the ferromagnet can lead 
to a second harmonic component in the CPR which is respon- 
sible for the finite supercuiTent at Tq^. The aim of the present 
work is to study the effects of the finite FS interface in the full 
range on the — tt crossover in diffusive SFS contacts. To 
our best knowledge there have not been studies of diffusive 
F Josephson junctions with the FS interface having an arbi- 
trary transparency. We show that even a small coiTection of 
the Josephson current due to weak disorder at the interfaces 
can produce a second-harmonic term which leads to a small 
finite supercurrent at the — tt as was observed in experiment. 

In this paper we investigate the effect of the disorder at 
the FS interface on the Josephson supercurrent in the dif- 



fusive F contact between two conventional superconductors. 
We adopt quasiclassical Green's function method in the dif- 
fusive limit implemented by the general boundary conditions 
of Nazarov^^, which allow us to obtain the Josephson current 
through the contact for an arbitrary strength of the barrier at 
the FS interfaces. We analyze the — tt transition and the 
CPR in terms of the exchange field and the thickness of the 
F contact in the full range of the FS interface transparency. 
In two limits of high- and low-transparent interfaces the — tt 
transition occurs on changing the temperature for different ex- 
change fields and thicknesses of the F layer For a transpar- 
ent interface the — tt transition does not occur for layers 
of thickness smaller than the coherence length £_ — -\/Co^imp 

[here ^o — ^'p /"'^olr — 0) is the superconducting coher- 

(s) 
ence length, Vp is the Fermi velocity in the superconductor, 

Ao (T = 0) is the amplitude of the superconducting order pa- 
rameter at zero temperature and i?imp is the mean free path in 
the F layer] even for very large exchange fields (throughout 
this paper we use the system of units in which the Plank and 
Boltzmann constants h ~ ks = 1). In contrast to this for a 
low-transparency interface the transition takes place for thin- 
ner contacts if the exchange field is strong enough. In these 
two limiting cases the current-phase relation is sinusoidal pro- 
vided that the weak- proximity approximation holds. This im- 
plies a zero supercurrent at the — tt transition. We show 
that the corrections to these results, both in high- and low- 
transparency cases, produce a second-harmonic term propor- 
tional to sin(2(/?). While for the low-transparency interfaces 
the second-harmonic term is so small at it can be neglected, 
it has some measurable value at the high-transparency limit 
which has oscillatory behavior as a function of the F thickness 
and the exchange field. For values of these quantities which 
lead to a positive value of the second-harmonic, the critical 
current at the — tt transition has a small finite value. This 
finding is consistent with the experiment'^, where evidence 
for the second harmonic was found for highly transparent in- 
terface samples. 

We further show that for interfaces far from being in the 
high- on low-transparency regime, the CPR deviates strongly 
from the standard sinusoidal form. The second-harmonic term 
at the — TT transition can be as the same order as the first 
harmonic for an intermediate value of the FS interface trans- 
parency, leading to a large residual supercurrent at the tran- 
sition, provided that the second-harmonic term is positive for 
the given exchange field and thickness of the contact. We an- 
alyze the finite supercurrent at the transition in the full range 
of the interface's transparency. 

In Sec. II, we introduce our model of a diffusive SFS con- 
tact. Assuming the weak-proximity effect regime we linearize 
the Usadel equations and the Nazarov boundary conditions at 
the FS interfaces with respect to the anomalous component of 
the quasiclassical Green's function in the ferromagnetic side 
of the interface. We present an expression for the Josephson 
current which is valid for an arbitrary barrier strength at the FS 
interfaces. Section III is devoted to analyzing the — vr tran- 
sition and the finite transition current in the full range of the 
parameters. We obtain analytical expressions for the second- 
harmonic corrections in the two limits of high and low trans- 
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FIG. 1 : Schematic of the SFS junction. The current flow between 
two superconductors (81,2) through a ferromagnetic layer (F). The 
FS interfaces have a finite transparency. 



parencies and also present the results for intermediate trans- 
parency . We also compare our results with the finding of the 
recent experiment. Finally we end with a conclusion in Sec. 
IV. 



II. THE JOSEPHSON CURRENT IN DIFFUSIVE SFS 
STRUCTURES 



In this section we calculate the Josephson current for a dif- 
fusive SFS junction in which the FS interfaces have an arbi- 
trary transparency. We consider the SFS junction schemati- 
cally shown in Fig. ^ It consists of two conventional super- 
conducting reservoirs Si and S2 which are connected by a thin 
ferromagnetic metal (F) layer The F layer has a thickness 
L and is characterized by a homogeneous spin-splitting ex- 
change field (energy) h . The Josephson current through the F 
layer is determined by the phase difference (p between the su- 
perconducting order parameters of the two superconductors. 
For a diffusive contact the thickness L is much larger than 
the mean free path of the impurity scattering, £i,„p. Band mis- 
matches and disorders at the interfaces lead to a backscattering 
of the quasiparticles. Thus the FSi and FS2 interfaces have 
finite transparencies which are determined by the interface re- 
sistances Rbg and R^j^ , respectively. We consider the limit of 
large exchange fields (h ^ Tc), when the superconducting 
correlations are strongly suppressed upon crossing from the 
superconducting side to the ferromagnetic side of the FS in- 
terfaces. In this case we can assume a weak superconducting 
proximity effect and apply a theory which is linearized with 
respect to the induced superconducting order parameter in the 
F layer. 

In the quasiclassical regime ii^p is much larger than the 
Fermi wavelength Ap which allows us to use the quasiclassical 
Green's function approach. For the diffusive junction we ap- 
ply the Usadel equation for the spin-resolved matrix Green's 

9<T fa 

n -9a 
flip scattering reads 



function c/a 



, which in the absence of spin- 



[{ujn + iah)Ty, + A(r),.9^] = Dd^g^drQa- 



(2.1) 



Here g„ and f„ are, respectively the normal and the anoma- 
lous components of the matrix Green's function for spin 



a {= ±1) electrons which depend on the Matsubara fre- 
quency ojn = TrT{2n + 1) and the coordinate r. The diffu- 
sion coefficient D — Vp ^imp/3, in which Vp is the Fermi 
velocity in the F layer and Ti (i=l,2,3) denote the Pauli matri- 

A(r) 
A*(r) 

is the superconducting pair potential (order parameter). The 
matrix Green's function g^ satisfies the normalization condi- 
tion: 



ces. The matrix A(r) 



, in which A(r) 



9l 



= 1, 



gl+un = i, 



(2.2) 



where f^{r, a;„, h) — /*(r, w„, —h). We have assumed a thin 
film structure with thickness L much smaller than the lateral 
dimensions, so that variations are only along the x axis in the 
direction perpendicular to the plane of the film. 

The current density is obtained from /^ through the relation 



Ji'P) 



iirT 
2epF 



E 



E 



=±1 w„ 



u 



dfl __ f^df, 



dx 



V 



dx 



(2.3) 



Here pp is the normal-state resistivity of the F layer. 

Inside the F layer A(x) = 0, but fa can still be nonzero 
due to the proximity effect. We neglect the suppression of the 
order parameter A close to the FS interface inside the super- 
conductors . Thus A(T) = Ao(r) exp (±i(^/2) is constant in 
the superconductors. With this assumption we solve Eq. ( 12. H 
inside the F layer and impose boundary conditions at the two 
interfaces which relates the values of the Green's functions at 
the superconducting and ferromagnetic sides of the two inter- 
faces. We apply the general boundary condition developed by 
Nazarov^^, which is valid for an arbitrary value of the barrier 
height at the interfaces. For g^ the boundary conditions are 
written as 



dga 
dx 



Rf 



M(0) 



(2.4) 



M{0) 



Rbo.L^ 

2To^Lmga(0-{L-)),ga{0+{L+))] 
4 + To,L ({5.(0-(L-)),5^(0+(L+))} - 2) ' 

(2.5) 



where Rp is the normal-state resistance of the F layer The 
transmission eigenvalues of the interfaces at 0, L are given by 
To.l(^) — cos^ 0/(cos^ ^ + 7b , ) with the angle 6 showing 
the direction of the Fermi velocity with respect to the axis nor- 
mal to the interfaces, "fb l — ^^b K^Pf) dimensionless 
quantities which measure the strength of the scattering from 
disorders at the interfaces 0, L, respectively, and A the area of 
the interfaces-^. The averaging over directions of the Fermi 
velocity is defined via 



M{e) 



j:i%d9 cosi0)M{9) 
j:i%d9 cos{e)T{e). 



(2.6) 



[ ] and {} show the commutator and the anticommu- 
tator in M{0), respectively. Also 9^(0^) [faiO^)] and 



.9(7 (O"*") [/cr(0"'")] are normal (anomalous) Green's functions 
at the interface in the superconductor and the ferromagnet 
sides, respectively. A similar notation holds for ga{L^) and 
fa-{L^) at the interface L. These boundary conditions are 
valid for arbitrary values of 76^ ^ . Equations ( I2.4t -( lz!6t are 
reduced to the boundary condition of KL in the limit of small 
transmissions (7b„ ^^ :§> 1) and to the condition of the continu- 
ity of the Green's function in the limit of high transparencies 

(76o,. -> 0). 

In a homogeneous bulk s-wave superconductor under the 
equilibrium condition, the solutions of the Usadel equations 
are found to be gs ~ ujn/^n, fs — A/Sl^, where ri„ = 
■\/Ag + 0^2.28 In the diffusive limit, the boundary values of 
the Green's functions in the superconductor sides of the inter- 
faces are approximated by the bulk values: 



5.(0-(L+)) = 9s = ^.. 



(2.7) 



fa{0-{L+)) = /,,, 






e^-^'-. (2.8) 



In general, the Usadel equation ( 12. U and the Nazarov bound- 
ary conditions Eqs. ( I2.4> -( |Z61 are nonlinear with respect to 
gcr and fa-. We linearize these equations by assuming a weak 
proximity effect in the F layer which implies \ga\ ~ 1 and 
\fa\ ^ 1. With this approximation Eqs. ( 12. Il l and ( 12. 5> are 
respectively reduced to the following relations 



(a;„ +iah)fa{x,ujn,h) - 



D d'^fa{x,UJn,h) 

2 dx^ 



= 0, (2.9) 



Af (61) = 



To,Lm9a{0-{L-)),ga{0+{L+))] 
2 + (gs - 1)To.l{9) 



± 



T^lWU (0-(i+)) [/.(0-(L-))/t(0+(L+)) 



[2 + (gs - l)To,Lm' 



(2.10) 



We note that although Eqs. ( I2.9l l and ( I2.10> are linear with 
respect to fa in the F layer, they are still nonlinear in Tq.l (d)- 
The anomalous Green's function fs in the superconducting 
side of the interface can be of order unity at low tempera- 
tures. Thus the FS interfaces have a very small width across 
which fa decreases from the bulk value fs to fp (fp is the 
anomalous Green's function in the ferromagnet side) which is 
assumed to be small at large exchange fields h ^ Tc- 
The solution of Eq. ( 12. 9> inside F layer has the form 



fa{x,UJn,h) ^ AnaCOSh{knaX)+ Bna'^Va^iknaX), (2.11) 



where kna = -J^iujn + icrh)/D. Replacing the solution Eq. 
(I2.11> in Eq. (12. 3> the expression for the Josephson current is 
obtained as: 

^(<^) = '-^77- E E KaIm{AnaBla)- 



epp 



(T — itl n— — OO 



(2.12) 



Here pp = l/(2e^A^o-D), where A'q is the density of states 
on the Fermi surface and Im{ ) denotes the imaginary part. 
The coefficients Ana and Bna are derived by applying the 



boundary conditions \2Al and ( I2.10> to the solution ( 12. IH and 
using Eqs. (12. 7> and (12. 8> . In this way we have obtained Ana 
and Bna and replaced the results in ( I2.12> to obtain the final 
expression for the Josephson current as 

n — + 00 

I{'P)^Iq^ X! Qna{^)si--a.tp, (2.13) 

a—±l n— — oo 

where Iq — 2TTTc/eRF, Tc is the superconducting critical 
temperature, and Qnaif) is given in the Appendix through 
Eqs. dAH - dXTt in which for simplicity we have assumed 

7&0 = 7bt = 76- 



0.0033 - 
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III. DISCUSSIONS AND RESULTS 

Equation ( I2.13> with Qna given by Eqs. (lAlt - JATt ex- 
presses the Josephson current through the diffusive F layer 
for an arbitrary value of 7b. In the limit of 7;, ^ 1 and 
7h ^ 1 the expression (I2.13> reduces to the results previously 
obtained for the cases of the fully transparent^ and the low- 
transparency'^ interfaces, respectively. In these two limits the 
current-phase relation is sinusoidal. Let us first find correc- 
tions to these results by expanding the Josephson current Eq. 
( I2.13> in terms of powers of 76 for small 7b and keeping the 
terms up to the order of 7;,. In the opposite limit 7^ ^ 1 we 
do expansion in the powers of 1/7^ and retain the terms up 
to the order of 1/7^ • The results in the limits of 76 ^ 1 and 



7h » 1 have, respectively, the forms 

I{ip) = /(I) sin if + /(2) sin(2(/5) 
/(V?) = J(i) sin(^ + J(2) sin(2(^). 



(3.1) 
(3.2) 



Here /('^•^J and j'^^^^^ are, respectively, the first- and the 
second-harmonic terms of the CPR in the limits of high and 
low transparency and are given in Eqs. ( IA8> -( lAT3t . 

Thus we find that a small deviation from the limits 76 = 
and the 7b ^ 1 has two effects on the Josephson current. 
First it contributes a small correction to the first-harmonic cur- 
rent term. Second it produces a small second-harmonic term. 
The second-harmonics terms Ina , Jna are oscillatory func- 
tions of h and L and depend on the temperature T. In spite of 
their smallness they can be important at the — tt transition 
where the first-harmonic terms vanish. While in the absence 
of the second-harmonic terms the critical current undergoes 

a continuous sign change at the crossover between and tt 

(2) (2) 
states, for nonvanishing and positive values of /„</ , Jna the 

sign change is discontinuous and the critical current jumps 
from a positive to a negative value. The discontinuous transi- 
tion is manifested as a nonzero minimum of the absolute val- 
ues of the critical current at the — tt transition temperature 

rp 25 

In the experimentally relevant limit when h ^ Tc and 
L > S^F (Cf — \/D/h is the superconducting correlation os- 
cillation length in the F layer) we find from Eqs. (IA8> -( IaTo1 
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FIG. 2: (a) Dependence of the Josephson critical current on the tem- 
perature for h/T^ = 9.1, L/C = 1.1, and 76 = 0.0016. Inset 
shows the temperature dependence near the transition temperature 
Tott/Tc = 0.134; /c(r) has a nonzero minimum at the transition 
temperature, (b) The current-phase relation for different tempera- 
tures around Tq.^. 



that always J^^^/J^^' < 10^^, so the second-harmonic term 
in the CPR may be neglected in the low-transparency limit. 
However for highly transparent interfaces Eqs. ( IA11> -( IaT31 
gives us that /(2)//(i) can be as large as lO^'^, a ratio that can 
be observed in experiment. 

Let us now compare our results with the findings of the 
experiment of Ref. 19. In this experiment the temperature 
dependence of the Josephson critical current was studied in 
the ferromagnetic alloy layer Cui_j;Nia; with x of the order 
0.48, connecting two Nb superconducting electrodes. The 
normal-resistance measurement of the junctions has shown 
that the samples have highly transparent FS interfaces. For 
two different thicknesses of the F layer Li = 17 nm and 
L2 = 19 nm, the — tt transition was observed at the tem- 
peratures 1.12 and 5.36 K, respectively. For the sample with 
L2 = 19 nm the critical current had a zero minimum, but 



they measured a nonzero supercurrent of about 4 /iA for the 
sample of Li — 17 nm at the — tt transition tempera- 
ture. The superconducting critical temperature of the samples 
Tc = 8.7 K and the exchange field was estimated to be of 
the order h — 5 — 10 meV. In Figs. |2]and|3]we present the 
Josephson critical current dependence on the temperature and 
the CPR of the corresponding two samples using our theoret- 
ical results given by Eqs. ( I2.13> and (IA8> -( lAT3t . We fix the 
critical temperature of the system, Tc, the Fermi velocity Vp of 
the superconductors, and the diffusion constant D of the two 
F layers with thicknesses Li 2 = 17 and 19 nm, by the values 
reported in the experiment'^. This leaves two fitting parame- 
ters to be obtained, which are the strength of disorder at the 
interfaces, 7^, and the exchange field h. To take into account 
the suppression of the order parameter at the superconduct- 
ing side of the FS interfaces, which is more pronounced for 
strong exchange fields h ^ Tc and a highly transparent inter- 
face, we have multiplied Aq{T) by a constant factor a = 0.5 
in our calculations. This, however, does not change the qual- 
itative behavior of /c with the temperature, but allows for a 
better fitting with the experimental curves. As in other theo- 
retical studies of the ferromagnetic Josephson junctions (see 
e.g., Ref. 9) the amplitude of our calculated critical current is 
substantially (~ 10'' times) larger than the measured curves. 
This can be attributed to the strong destruction of the Andreev 
pairs due to the spin-flip scattering induced by the superpara- 
magnetic behavior of the Cu-Ni alloy"', which is not taken 
into account in our calculation. 

For the thickness 17 nm the results are presented in Figs. 
Ha) and 2(b). We have taken D = 5 cmVs and h/Tc = 9.1. 
In Fig. |3a) the temperature dependence of the critical cur- 
rent is presented for one value of 7^. There is a good agree- 
ment with the experimental curve at 7^ = 0.0016. Both the 
functional dependence and the transition temperature are re- 
produced correctly. At the transition temperature To^/Tc = 
0.134 the critical current has a finite value Iq/Iq ~ 3 x 10~^. 
The finite critical current results from the appearance of the 
second-harmonic Josephson current in the CPR. This is shown 
in Fig. |2jb) in which we present the CPR for different temper- 
atures around the transition when 7b — 0.0016. Approaching 
the transition from lower and higher temperatures the CPR de- 
viates from a sinusoidal form and tends to be a purely second 
harmonic proportional to iim{2ip) at Tott- 

Figure, ^a) shows the temperature dependence of the crit- 
ical current for the second thickness 19 nm. The diffusion 
constant and the exchange field have the values D ~ 5 cm^/s 
and h/Tc = 7.93. The value of 7^ = 0.002 corresponds to the 
transition observed in the experiment with Tqt^/Tc — 0.619 
and a zero current minimum. The temperature dependence of 
the CPR is presented in Fig. |3jb) for 7^ = 0.002. In con- 
trast to the Li = 17 nm case the CPR remains almost sinu- 
soidal crossing the transition temperature. At the transition 
temperature I{(p) = for all the phases, which represents a 
continuous — tt transition. 

Let us finally analyze the — tt transition at intermediate 
transparencies. In Fig. ^a) the temperature dependence of 
the critical current Ic{T) is presented for different 7^ when 
h/Tc = 5 and L/£^ = 1.1. For these values of h and L the 
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FIG. 3: The same as Fig. Shut for h/Tc = 7.93, L/^ = 1.23, and 
7b = 0.002. 



— TT transition does not take place for fully transparent inter- 
faces, 7b = 0. At a lower transparency 7b = 0.52 the transi- 
tion occurs at Tqt^/Tc = 0.113 with a large residual minimum 
supercurrent. On decreasing the transparency, Tott increases 
toward the critical temperature Tc and the transition current 
decreases and vanishes as Tott ^ Tc- Thus we find that the 
highest values of the second-harmonic term can be achieved 
at the intermediate transparency. This can also be seen from 
the CPR which is presented in Fig. ^h) for 7b = 0.52 at dif- 
ferent temperatures. There is a strong deviation from a sinu- 
soidal CPR relation and I{lp) oc sin(2(p) exactly at the tran- 
sition temperature. We also obtained that for a fixed thick- 
ness by increasing the exchange field the transition tempera- 
ture Tott shifts toward Tc and the minimum current amplitude 
decreases. Similar behavior holds on increasing the thickness 
and fixing the exchange field. As a result the — tt transition 
taking place at a higher temperature has a smaller minimum 
current. We expect that even for the sample of i = 19 nm 
there is an extremely small minimum supercurrent which was 
not observable in experiment. 
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FIG. 4: (a) The Josephson critical current versus T/Tc when h/Tc — 
5 and L/^ = 1.1, for different 7;,, in the intermediate transparency 
regime, (b) The current-phase relation dependence on the tempera- 
ture for 7i, — 0.52. 



IV. CONCLUSION 

In conclusion we have explained that the recently observed 
finite supercurrent at the — tt transition temperature of diffu- 
sive feiTomagnetic Josephson junctions can originate from a 
weak quasiparticle backscattering due to small imperfections 
at the FS interfaces. Using the quasiclassical Green's func- 
tions approach in the diffusive limit implemented by the gen- 
eral boundary conditions of Nazarov at the FS interfaces we 
have obtained an expression of the Josephson current through 
the F layer which is valid for any strength of the interfacial 
barrier. At the limiting cases of a perfectly transparent and a 
low-transparency interface our results reduce to the previously 
obtained expressions of the Josephson currents with sinu- 
soidal cuiTent-phase relations. We have found that in both lim- 
its the coiTections to these results contain a second-harmonic 
component. While for the experimentally relevant values of 
the ferromagnetic layer thickness {L/^p 3> 1) and the ex- 
change field {h/Tc ^ 1) the second-harmonic correction is 



very small and can be neglected in the low-transparency limit, 
it can have an observably large positive value for highly trans- 
parent interfaces. We have demonstrated that in the highly 
transparent interfaces the positive second-harmonic current 
results in a finite supercuiTent at the — tt transition, which 
has been observed in experiment. We have compared our re- 
sults with the experimental curves and found a satisfactory 
agreement. We also have shown that for an intermediate trans- 
parency of the interfaces the second-harmonic Josephson cur- 
rent can be as large as the first harmonic, leading to a large 
residual supercuiTent at the — tt crossover in this limit. 
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APPENDIX A: THE EXPRESSION Q„<, IN THE 
JOSEPHSON CURRENT FOR THE CASE -yt^ = jb^ = 7i> 

After deriving the coefficients Ana and Bna- and replacing 
the results in Eqs. ( I2.12> we obtain the final expression for 
the Josephson cuiTent which is given by Eqs. ( 12.131 1. in which 
Qnai^p) is given hy 



Pi 
P2 



Ps = 



32b''LT k^apiAl 

Tc Po 

x[/3i + 16cA2(-7,^cosv? + /32)]', (A1) 
{Pl + 32cAl{pip2 
+AcAl[-Ka7b^^ cos(2(p) + P3]})\ (A2) 

fi^ [iQbgsknalbS, coah{kncrL) 

+ (6462^2 ^ fc2^^2^2) sinh(fc„,i)] , (A3) 
kna^b^cosh{knaL) + Sfo^s sinh(fc„c,L), 

(A4) 
(646^52 ^ fc2^^2^2) cosh(2fc„,i) 

kncrlbi. 



+6Ab^gU-l 

no) 



Ahg, 



'i9s-i)ne) 
T\e) 

\{gs-l)T{e) + 2) 



sinh(2fc„^L),(A5) 
(A6) 
2). (A7) 



In the limits of 76 <C 1 and jb ^ 1, we can expand the 
Josephson current Eq. (I2.13> in terms of jb- In the limit of 
7b <C 1, we keep terms up to the order of 7^. In the opposite 
limit of 7f, ^ 1 we do the expansion in powers of l/7f, and 
retain terms up to the order of 1/7^. As a result, the Josephson 
current can be expressed in terms of the first and second har- 
monics whose coefficients in the limits of 7^ <C 1 and 75 ^ 1 



are /f^^^) ^jj^j j(i,2)^ respectively, 






j(2) 



n^+c>o 



E E 

(T— ±1 n — — CX2 



T,w2p2 sinh(fc„,L) ' 



Fi = l--/b^RiAf,cothikn,L) 



P2 = 
Ri = 



A2 

^0 



ids +9^1 ' 



(A8) 



i?2 



i?2 = 3 



224(1 + ,g,)' [1 + 3w„a.(l + 9s)/i7Al)] 



^l^l 



n— +OC 



E E 



(T— ±1 n— — 00 



-27bg£rfc,2,,A4 



(A9) 
(AlO) 



i?4 

J(2^ 



n— +00 



LTAIR3 



_ Rj coth(fc„cr-^) 



(All) 



2 , A2[(l + g,)-A2] 2g,(A:„,0"' 



(1 + .gs)2 sinh(A;„cri) ' 

(A12) 



n— +00 



-LTA^ 



.til nf^'oo 4(1 + gs)Tejle^tkl, Smh\kn.L) ' 

(A13) 
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